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ABSTRACT 

Starburst galaxies are primary feedback sources of mechanical energy and metals, which 
are generally measured from associated X-ray emission lines assuming that they are from 
the thermal emission of the outflowing hot gas. Such line emission, however, can also arise 
from the charge exchange X-ray emission (CXE) between highly ionized ions and neutral 
species. To understand the feedback of energy and metals, it is crucial to determine the origin 
of the X-ray emission lines and to distinguish the contributions from the CXE and the thermal 
emission. The origin of the lines can be diagnosed by the Ka triplets of He-like ions, because 
the CXE favors the inter-combination and forbidden lines, while the thermal emission favors 
the resonance line. We analyze the triplets of O vn, Neix, and Mgxi observed in the XMM- 
Newton reflection grating spectra of the starburst galaxy M82. The flux contribution of the 
CXE is 90%, 50%, and 30% to the O vn, Ne ix, and Mg xi triplet, respectively. Averaged over 
all the three triplets, the contribution of the CXE is ~ 50% of the total observed triplet flux. To 
correctly understand the hot outflow of starburst galaxies, it is necessary to include the CXE. 
Based on the measured CXE contributions to the O vn, Ne ix, and Mg xi triplets, we estimate 
the relative abundances of O, Ne, and Mg of the outflow and find they are similar to the solar 
ratios. 
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1 INTRODUCTION 

Starburst galaxies are the primary sources that eject energy 
and metals into the intergalactic medium. The properties of the out- 
flow are fundamental to the understanding of the feedback process 
of galaxies. The energy and metals of the outflow are generally 
measured from the X-ray emission, particularly emission lines, as- 
sociated with starburst galaxies. In past studies, the X-ray emission 
lines are assumed to come from the hot outflow itself, the tempera- 
ture and metal abundances of which are extracted based on thermal 
models. 

The X-ray line emission, however, can arise not only from 
the hot gas, but also from the interaction between the hot gas and 
the neutral cool gas. For example, when the solar wind interacts 
with a comet, the highly ionized ions in the wind can readily pick 
up electrons from the neutral species of the comet. The product 
ions remain highly ionized and are left in excited states. They will 
emit X-ray photons when they decay to ground states. This pro- 
cess is called charge exchange X-ray emission (CXE, also called 
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charge transfer) and explains the bright cometary X-ray emission 
(e.g. Lisse et al. 1996; Cravens 1997, 2002). It can be represented 
by: 



A" + + N -> A u '-' )+ ' + N + , 



(1) 



where a highly ionized ion A' l+ (like O Vin, Ne x) picks up an elec- 
tron from a neutral species N (like H, H2), producing an exited 
ion A <9_1)+ *, which will emit X-ray photons when it decays to the 
ground state. For the historical studies and recent developments of 
the CXE, we refer to Dennerl (2010) and references therein. 

Different from the thermal emission, the CXE contributes only 
emission lines. If the X-ray line emission, or part of it, is due to the 
CXE, the measurement of the thermal and chemical properties of 
the hot outflow based on thermal-only models will be misleading. 
To correctly understand the hot outflow, it is crucial to reveal the 
origin of the X-ray emission lines and to distinguish the contribu- 
tions from the CXE and the thermal emission. 

Based on the tight spatial correlation of Ho- and X-ray emis- 
sions, Lallement (2004) speculated the importance of the CXE be- 
tween the hot outflow and the cool halo gas for starburst galaxies. 
From Suzaku and XMM-Newton observations of the cap above the 
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disk of M82, Tsuru et al. (2007) showed a marginal detection of an 
emission line at 0.459 kev, which may be due to the CXE of CVI 
(n = 4 — > 1). From XMM-Newton EPIC spectra of the central re- 
gion of M82, Ranalli et al. (2008) also reported two lines (around 
10 and 16 A), which can not be adequately accounted for by ther- 
mal models and may be contributed by the CXE from neutral Mg 
and O vm. 

One diagnostic that can be used to determine the origin of the 
X-ray emission lines is the Ka triplet between n = 2 shell and n = 1 
ground state of He-like ions, which consist of a resonance line, two 
inter-combination lines, and a forbidden line (for a recent review, 
see Porquet et al. 2010). For a thermal plasma, the electron col- 
lisional excitation is efficient and favors the resonance line, while 
for the CXE, the de-excitation favors the triplet states and thus the 
inter-combination and forbidden lines. Therefore, the line ratios of 
the triple can determine the origin of the X-ray emission lines. If the 
forbidden line dominates, the CXE will be the most feasible origin, 
because the hot plasma is already highly ionized and photonization 
is negligible. 

In this letter we analyze the Ka triplets of He-like ions of 
O vn, Ne IX, and Mg XI of M82 observed by XMM-Newton Reflec- 
tion Grating Spectrometers (RGS, den Herder et al. 2001). M82 
is one of the nearby brightest galaxies in X-ray and regarged as a 
prototypical starburst. A bipolar outflow originated from the nu- 
clear starburst extends beyond 5' (about 5 kpc) from the disk (e.g., 
Stevens et al. 2003). Ranalli et al. (2008) have analyzed the XMM- 
Newton RGS O vil triplet of M82 and they fitted three Gaussians 
to the blended triplet. They showed that the line intensity ratios are 
marginally consistent with the CXE. In this letter, we decompose 
the triplets into the CXE and the thermal emission using the labo- 
ratory measured CXE spectrum of O vil triplet. As a result, we can 
quantify the contribution of the CXE to the observed Ka triplets 
of He-like ions, including O vn, Ne IX, and Mg XI. The RGS also 
allow us to study the O vil triplet at different spatial regions. 

The local emissivity of the CXE of one particular line (assum- 
ing single electron capture) P c oc n A q+ Vcrn neu , where n A q+, n„ eu , V, 
and o~ denote the ion density, the neutral density, the relative ve- 
locity, and the charge exchange cross section, respectively. If we 
know the CXE fluxes of the O VII, Ne IX, and Mg XI triplets, we can 
estimate the relative abundance of O, Ne, and Mg with the temper- 
ature information. It provides important insight to the composition 
of the hot outflow, especially when the assumption of the thermal- 
only origin of the X-ray line emission is poor. We will make a first 
attempt to estimate the relative abundances of O, Ne, and Mg using 
the CXE fluxes of their He-like triplets. 

The paper is structured as follows. We describe the observa- 
tion data in §2. The analysis results are presented in §3. Our con- 
clusion and discussion are given in §4. Errors are given at 1 cr con- 
fidence level. 



2 OBSERVATION DATA 

We use four archival XMM-Newton RGS observations of M82, 
the obs-IDs of which are 0112290201, 0206080101, 0560590201, 
and 0560590301. The total effective exposure time is ~ 120 ks af- 
ter removing intense flare periods, which doubles the useful dataset 
used by Ranalli et al. (2008). The most recent version of Science 
Analysis System (SAS, 10.0) is used for the reduction of photon 
events. In these observations, the RGS slit is generally not parallel 
to the outflow direction and the extracted emission lines are broad- 
ened due to the tilt of the outflow. Thus we limit our extraction to 
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Figure 1. Fluxed XMM-Newton RGS spectrum of M82 extracted from the 
central 30" radius of the cross-dispersion direction. Both RGS1 and RGS2 
are combined together. The spectrum is dominated by Lya lines of H-like 
and He-like states of O, Ne, and Mg and neon-like Fe lines. The triplets of 
O VII, Ne IX, and Mg XI are analyzed below. 



events within 30" radius of the cross-dispersion direction. We also 
extract the Ovil triplet from -90" to -30" and from 30" to 60" 
of the cross-dispersion direction, where the O vil triplet is bright 
enough for further study. For the region larger than 60", the O vil 
triplet is heavily blended by the structures of the X-ray emission 
and does not allow a simple fitting as we do below, so it is excluded. 

There are two RGS on the XMM-Newton . CCD4 in RGS 2 
covering the O vil triplet failed, and CCD7 in RGS 1 covering the 
Ne IX triplet also failed. Thus, the O vil triplet data are from RGS 1 
and the Ne IX triplet data from RGS2, while the Mg XI data are from 
both RGS1 and RGS 2. The Sixm triplet around 6.7 A is neglected 
here as the effective area is small and the photon counts are too few 
to obtain robust results. 



3 RESULTS 

The observed XMM-Newton RGS spectrum of M82, corrected 
for the effective area, is shown in Figure 1. As can be seen, the 
emission is dominated by Lya lines of H-like and He-like states 
of O, Ne, and Mg and neon-like Fe lines. Below we analyze the 
Ka triplets of He-like ions of O vn, Ne IX, and Mg XI. The intensity 
ratio of resonance, inter-combination, and forbidden lines can tell 
the origin of the X-ray emission lines. 



3.1 O vil triplet 

Let's first focus on the observed O vil triplet within the cen- 
tral -30" - 30" region, which is plotted in the left panel of Figure 
2. A continuum interpolated between the line wings has been sub- 
tracted. It shows clearly that the forbidden line is stronger than the 
resonance line. The dominance of the forbidden line can not result 
from the thermal excitation, but is a signature of the CXE. The lab- 
oratory measured CXE of O vn triplet using the spare X-ray micro- 
calorimeter from ASTRO-E (Suzaku) mission (Beiersdorfer et al. 
2003) is over-plotted as dashed line in Figure 2. The laboratory 
measured CXE of O vn triplet fits the observed O vil triplet well. 

To quantify the contributions of the CXE and the thermal 
emission, we fit the observed triplet with two components of the 
CXE and the thermal emission. The CXE component consists of 



Charge Exchange X-ray Emission ofM82 3 



2.0x10"* 




21-0 21.5 22.0 22.5 21.0 21.5 22.0 22.5 21.0 21.5 22.0 22.5 

A (A) A (A) A (A) 



Figure 2. RGS spectra of the O VII triplet in three extraction regions. A continuum, which is calculated as the interpolation between the line wings, has been 
subtracted. Laboratory data for the CXE of the O VII triplet by Beiersdorfer et al. (2003) is over-plotted as the dashed line in the left panel. The RGS spectrum 
agrees well with the CXE data. The solid line is the sum of the fitted model of the CXE (dotted line) and the thermal emission (dot-dashed lines). Most of the 
flux of O VII triplet is due to the CXE. 



Table 1. Results of fitting the CXE and thermal emission to the triplets of O VII, Mg XI, and Ne IX. 
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Ovn(-30" - 


30") 


8.6+1.3 


1.0+1.2 


0.12±0.02 


0.14±0.01 


7.0/6 


21.6 


21.8 


22.1 


0.624 


0.070 


0.306 


OVII(-90" - 


-30") 


4.9±0.9 


0.1 ±0.9 


0.15±0.01 


0.11±0.01 


6.4/6 


21.6 


21.8 


22.1 


0.624 


0.070 


0.306 


OVII(30" - 


60") 


2.3±0.8 


0.3±0.8 


-0.1 ±0.05 


0.15±0.03 


5.2/6 


21.6 


21.8 


22.1 


0.624 


0.070 


0.306 


MgXI 




2.6± 1.4 


7.1+1.5 


0.03±0.006 


0.05±0.005 


11.1/9 


9.17 


9.23 


9.31 


0.587 


0.101 


0.312 


NeIX(f = 


1.5*) 


5.1+1.2 


5.7±1.0 C 


0.06±0.007 


().06±0.006 


10.2/8 


13.45 


13.55 


13.7 


0.138 


0.027 


0.076 


NeK(g = 


4*) 


4.6±1.4 


6.5±1.2 C 


0.06±0.008 


0.06±0.()07 


11.5/8 


13.45 


13.55 


13.7 


0.262 


0.052 


0.144 



Note: for the meaning of symbols, see eq. (2); "f c and /, are in units of photons/s/cm 2 ; *^=f is in units of the solar ratio; c the value /, of NelX is for the 
thermal emission of both Ne IX triplet and Fe XIX lines. 



three Gaussians centered on the triplet lines, the ratios of which 
are fitted to the laboratory measured CXE triplet by Beiersdorfer 
et al. (2003). The thermal component also consists of three Gaus- 
sians, the ratios of which are calculated at T = 8 x 10 6 K, which 
is the estimated temperature of the hot gas (see §3.2). The modeled 
spectrum can be written as 



fill: 



1 



(A - Xj - AA) 



2a> 



(2) 



where Aj is the wavelength of resonance, inter-combination, and 
forbidden lines, c ; the corresponding ratios for the CXE, r, the 
ratios for the thermal emission, f c and /, is the flux of the CXE 
and the thermal emission respectively. The CXE ratios c rj ^=0.317, 
0.149, 0.535; the thermal ratios t riJ are calculated using APEC 2.0 
(Smith et al. 2001) and listed in Table 1. We assume that the disper- 
sion cr A and the wavelength shift AA are the same for both compo- 
nents. Since the observed RGS lines are dominated by the spatial 
extent of the source, it means that we assume both components are 
from a similar spatial volume. The fitted results are listed in Table 
1 and plotted in Figure 2. The thermal component only accounts 
for 10±12% of the observed Ovil triplet. As the CXE dominates, 
the results are insensitive to the adopted temperature. For example, 
if we adopte T = 5 x 10 6 K, the thermal fraction is 11±14%. For 
higher temperatures, the thermal ratios t T jj are almost unchanged 



and the results are not affected. Thus we conclude that most of the 
observed O VII triplet is due to the CXE. 

We can estimate the CXE contribution to the observed O vm 
Lya- (19 A) flux within the same region, which is ~ 1.8 x 10~ 4 
photons/s/cm 2 . In §3.2, we will see that the flux ratio of Mgxn 
Lya to Mgxi triplet indicates a temperature around 8 x 10 6 K. At 
such a high temperature, the ion fraction of O vm is negligible and 
O IX ions dominate charge exchange collisions. While it is possible 
for an O IX ion to capture two electrons to become an O VII ion di- 
rectly, the cross section of double-electron charge exchange is gen- 
erally lower than 20% of that of single-electron charge exchange 
(Greenwood et al. 2001). Thus we consider the case of multiple 
charge exchange collisions. That is, an Oix ion experiences one 
charge exchange collision becoming an Ovm ion, then becomes 
an Ovil ion through another collision. As the transition ratio of 
n(Ss 3 — » 1) to n(= 2 — > 1) is measured to be similar for the col- 
lision of Oix and Ovm (Greenwood et al. 2001), in the case of 
multiple collisions, the emission of an O vil-triplet photon is cou- 
pled with the emission of an O vm-Lya photon. In other words, the 
CXE flux of O vm Lya is the same as that of O vn triplet. Since 
the fitted CXE flux of O vil triplet is 8.6 x 10~ 5 photons/s/cm 2 , we 
expect that about half of the total observed flux of O vm Lya is due 
to the CXE. The CXE contribution to the O vm Lya emission will 
be larger if the assumption of multiple collisions is not fulfilled. 
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Figure 3. RGS spectrum of the Mg XI triplet. An interpolated continuum 
has been subtracted. The solid line is the sum of the CXE component (dotted 
line) and the thermal component (dashed line). The Mg XII Lya line is also 
plotted. 



Figure 4. RGS spectrum of the Ne IX triplet. The continuum has been sub- 
tracted. The NelX triplet is mixed with FeXIX lines. The solid line is the 
sum of the CXE (dotted line), the thermal Ne IX triplet (dashed line), and 
the thermal FeXIX emission (dot-dashed line). 



The observed O VII triplets from -90" to -30" and from 30" 
to 60" are plotted in the middle and right panels respectively. They 
also show the dominance of the forbidden line. That is, the CXE 
occurs not only at the starburst nucleus (~ 30"), but also at the 
outflow regions. The fitted results are listed in Table 1. The different 
A/1 reflects the tilt of the outflow relative to the cross-dispersion 
direction mentioned in §2. The observed Ovm Lya fluxes of the 
two regions are 8.6 and 4.9 x 10~ 5 photons/s/cm 2 respectively. The 
observed flux ratios of O VII triplet to O vm Lya are similar for all 
three regions, which implies a same scenario of multiple charge 
exchange collisions. 



3.2 Mg XI triplet 

The observed Mgxi triplet and Mgxn Lya line are plotted 
in Figure 3. Different from the Ovu triplet, the Mgxi triplet is 
dominated by the resonance line, i.e., by the thermal emission. 
This enables an estimation of the temperature of the hot gas as- 
suming there is no CXE. The ratio (0.57) of the observed flux of 
Mgxn Lya line (5.6 x 10~ 5 photons/s/cm 2 ) to that of Mgxi triplet 
(9.7 x 10~ 5 photons/s/cm 2 ) implies a temperature of 8 x 10 6 K. We 
fit the observed Mg xi triplet with both the CXE and the thermal 
emission. The ratios c r>1> / of the CXE are assumed to be the same 
as those of the laboratory measured CXE of the O VII triplet. The 
ratios of the thermal component are calculated at T = 8 x 10 6 K. 
The fitted results are plotted in Figure 3 and listed in Table 1 . The 
contribution of CXE is about 30± 15% to the observed Mg xi triplet. 

We can perform a self-consistency check on the temperature 
estimation by taking into account the CXE. The ion fraction of 
Mgxm and Mgxn at T = 8 x 10 6 K is 0.28 and 0.44 respectively. 
In the case of multiple charge exchange collisions, the CXE flux 
of Mg xil Lya line is 40% of that of Mg xi triplet. Excluding the 
CXE flux, the remaining thermal flux of Mgxn Lya is 4.6 x 10~ 5 
photons/s/cm 2 . The ratio of the thermal flux of Mg xil Lyo- to that of 
Mgxi triplet (7.1 x 10~ 5 photons/s/cm 2 ) is 0.65, which corresponds 
to a temperature of ~ 8.2 x 10 6 K. If single collision happens, the 
thermal flux of Mgxn Lya is 4 x 10~ 5 photons/s/cm 2 , correspond- 
ing to a temperature of ~ 8x 10 6 K. Thus the temperature estimated 
without the CXE is consistent with that considering the CXE. 



3.3 Ne ix triplet 

The observed Ne IX triplet is plotted in Figure 4. It is similar 
to the Mgxi triplet (right part of Fig. 3). The NelX triplet, how- 
ever, is mixed with Fexix lines. At T = 8 x 10 6 K, the Fexix 
lines are comparable to, or even stronger than, the Ne IX triplet. We 
take into account the Fe xix lines by including several other Gaus- 
sians centered on strong Fe xix lines. The abundance ratio of Ne/Fe, 
obtained by Ranalli et al. (2008) from fitting the RGS spectrum as- 
suming thermal-only models, is 1 .5 times the solar value (calibrated 
to Lodders 2003). However, their fitted Mg/Fe ratio is 4 times the 
solar value. Since Fe is mainly produced in SNe la and a-elements 
are enriched by SNe II, in principle, a/Fe should be similar. Thus 
we adopt two values (1.5 and 4) for the Ne/Fe ratio. The fitted re- 
sults are listed in Table 1 and the result of the Ne/Fe ratio of 1.5 is 
plotted in Figure 4. 

We see from Table 1 that the fluxes of the CXE do not change 
much for the Ne/Fe ratios of 1.5 and 4. This is because the mixing 
of Fe xix lines is mainly around the resonance line, the Ne/Fe ratio 
does have much effect on the relative importance of the resonance 
and forbidden lines. The CXE flux of the NelX triplet (~ 5 x 10~ 5 
photons/s/cm 2 ) dominates over the thermal flux of the Ne IX triplet 
(1.4 and 3xl0~ 5 photons/s/cm 2 for Ne/Fe=1.5 and 4 respectively) 
and accounts for ~ 50% of the total observed flux. 



4 CONCLUSION AND DISCUSSION 

We have analyzed the Ka triplets of O VII, Ne IX, and Mg XI 
of the XMM-Newton RGS spectra of M82. We show that most of 
the observed O vil triplet is due to the CXE. The contribution of the 
CXE is 30% to the Mg xi triplet. Taking into account the mixing of 
Fe xix lines, the fitted CXE of Ne IX triplet dominates over its ther- 
mal emission. Averaged over all the three triplets, the contribution 
of the CXE is about 50% to the total observed flux. 

As stated in the introduction, the local emissivity of the CXE 
of one particular line P c = (£ ? f q n Aq +)Van neu , where the sum is 
over all possible ionization states that can charge exchange to pro- 
duce the line with an efficiency f q . We can estimate the relative 
abundances of O, Ne, and Mg based on the CXE fluxes of their Ka 
triplets. We consider the case of multiple charge exchange colli- 
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sions as discussed in §3.1. An additional complication is the fore- 
ground absorption. While the column density N H can be obtained 
from fitting the whole spectrum as discussed below, here we adopt 
a column density N H = 0.2x 10 22 cnT 2 , which is taken from Ranalli 
et al. (2008). Then the absorption factor (/„;,.,) relative to the O VII 
triplet is 3 and 4.5 for the Ne IX and Mg XII triplet respectively. 

We use the relation of the charge exchange cross section 
<x cc Z, where Z is the element number (e.g. Wargelin et al. 2008). 
We assume the relative velocity is the same for all three elements as 
their element numbers are similar. We also assume the efficiency of 
producing the triplet f q is the same for all three elements. The ob- 
served CXE flux of the triplet f c cc n A f A Zf abs , where n A is the num- 
ber density for element A and f A the fraction of both fully stripped 
and H-like ions for element A. This leads to n A oc f c /(fAZf a bs)- 
Using the fitted f c in Table 1 , we find that the chemical abundance 
relative to O is 0.16 and 0.06 for Ne and Mg respectively. These ra- 
tios are similar to the solar ratios of 0.15 and 0.07, calculated from 
the solar abundance table of Lodders (2003). 

The O abundance of the hot gas of M82 is found to be about 4 
times smaller than other a elements in thermal-only models (Read 
& Stevens 2002; Origlia et al. 2004), which is difficult to explain. 
However, the abundance ratios of O, Ne, and Mg estimated from the 
CXE are similar to the solar ratios. Therefore, the low abundance of 
O obtained from thermal-only models is likely to be a consequence 
of neglecting the CXE. 

As the average contribution of the CXE to the triplets is ~ 
50% and the contribution differs for different lines, the CXE will 
affect the estimation of the thermal and chemical properties of the 
hot gas. In previous studies of M82 (e.g., Read & Stevens 2002; 
Ranalli et al. 2008), multiple-temperature models have been used 
to fit the observed spectra. The multiple temperatures are likely due 
to the neglect of the CXE. For example, if neglecting the CXE, the 
temperature estimated from the observed flux ratio of O vm Lya to 
O vil triplet will be 3.5 x 10 6 K, which is misleading. To properly 
estimate the temperature and metal abundances, it is necessary to fit 
the whole observed spectrum with both the CXE and the thermal 
emission. For such a task, a detailed model of the spectra of the 
CXE, including the CXE of Fe xvn is needed. We plan to return to 
this in a later work. 
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